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\ Abstract. A general type of localized excitations, folded solitary waves and foldons, are 

defined and studied both analytically and graphically. The folded solitary waves and foldons 
may be "folded" in quite complicated ways and possess quite rich structures and abundant 
interaction properties. The folded phenomenon is quite universal in the real natural world. 
The folded solitary waves and foldons are obtained from a quite universal formula and the 
■ universal formula is valid for some quite universal (2+l)-dimensional physical models. The 

"universal" formula is also extended to a more general form with many more independent 
arbitrary functions. 
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G^ . In the study of nonlinear science, soliton theory plays a very important role and has 

O ■ 

O ' been applied in almost all the natural sciences especially in all the physical branches such 

o : 

as the fluid physics, condense matter physics, biophysics, plasma physics, nonlinear optics, 
(N ■ 

O ' quantum field theory and particle physics etc.. Almost all the previous studies of soliton 

theory especially in high dimensions are restricted in single valued situations. However, 
Ch ■ the real natural phenomena are very complicated. In various cases, it is even impossible 

to describe the natural phenomena by single valued functions. For instance, in the real 
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■ natural world, there exist very complicated folded phenomena such as the folded protein 

folded brain and skin surfaces and many many other kinds of folded biologic systems [0. 
The simplest multi-valued (folded) waves may be the bubbles on (or under) a fluid surface. 
Various kinds of ocean waves are really folded waves also. 

To study the complicated folded natural phenomena is very difficult. Similar to the single 
valued cases, the first important question we should ask is: Are there any stable multi- 
valued (folded) localized excitations? For convenience later, we define the multi-valued 
localized excitations folded solitary waves. Furthermore, if the interactions among the folded 
solitary waves are completely elastic, we call them foldons. In (l+l)-dimensional case, the 
simplest foldons are so-called loop solitonsj^ which can be found in many (l+l)-dimensional 

integrable systems Q and have been applied in some possible physical fields like the string 
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interaction with external field[0, quantum field theory and particle physics0. However, 
in our knowledge, there is no study at all for the possible higher dimensional foldons. That 
is the main topic of this letter. 

In (l+l)-dimensions, it has been proven that when a physical model can be expressed 
by a partial differential equation (PDE), then under some suitable approximations, one 
can always find nonlinear Schordinger type equations^ that is why the (l+l)-dimensional 
soliton theory can be successively applied in almost all the physical branches. Similarly, if 
(2+l)-dimensional physical model can be expressed by a PDE, then under some suitable 
approximations, one can find the Davey-Stewartson (DS) type equations]^. 

Recently, it is found that a quite "universal" formula P 

U = - ■ ■ ■ r^, A = aoa3-aia2, (1) 

(oo + aip + a2q + aspq)^ 

is valid for suitable fields or potentials of various (2+l)-dimensional physically interesting 
integrable models including the DS equation, the dispersive long wave equation (DLWE)||10||, 
the Broer-Kaup-Kupershmidt (BKK), and so on. In Eq. (|l]), p = p{x, t) is an arbitrary 
function of {x, t}, q = q{y, t) may be either an arbitrary function for some kinds of models 
such as the DS system or an arbitrary solution of a Riccati equation for some others, say, the 
DLWE while oq, cii, 02 and 03 are taken as constants. One of the most important results 
obtained from (|l|) is that for all the models mentioned above there are quite rich localized ex- 
citations such as the solitoffs, dromions, lumps, breathers, instantons, ring solitons, peakons, 
compactons, localized chaotic and fractal patterns and so onp[. 

Now the natural and important question is: Can we find some types of (2+l)-dimensional 
foldons from the "universal" formula? Fortunately, the answer is obviously positive because 
of the arbitrariness of the functions p and q included in the universal formula. To find more 
general types of foldons, we extend the universal formula to a more general form for the 
DLWE before to give out concrete foldons. For the (2+l)-dimensional DLWE 

Uyt + Vxx + UxUy + UUxy = 0, Tjt + + TjU^ + Urjn, + Un,xy = 0, (2) 

(|I|) is valid for the field v = —rj — 1 and the related exact solution reads 

1 rr , 2p^{ai + a3q) 

V = —1] — 1 = U, u = ± h Mo, (3) 

ao + aip + a2q + a^pq 

with p being an arbitrary functions of {x, t}, q = q{y,t) being an arbitrary solution of 

the Riccati equation, qt — aoCo — (aiCi + a2Co — aoC2)q — (oaCi — 0202)5^ = 0, and uq = 
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~Px^[Pt^Pxx — doCi — (fliCi + a2Co + aoC2)p — (aiC2 + asCo)^^]. The (H-l)-dimensional DLWE 
{y = X of (^) is also called the classical Boussinesq equation. There exist a large number 
of papers to discuss the possible applications and exact solutions of the (l+l)-dimensional 
DLWE[rT|. Various interesting properties of the (2+l)-dimensional DLWE have been studied 



by many authors|T^-[|T^. 

To extend the universal formula to a more general form, we use the multi-linear variable 
separation approach again. For the DLWE (0), the transformation {u = Uq ± 2^, rj = 
2^ — — 1} with Mo = uo(a;,t) degenerates two equations of (H) to a trilinear form 

[fxxxy i (fxyt + UQxfxy + Wq/xxj/)]/^ ~ [{fxxx i Mq/xx i fxt)fy i fxfyUQx 

+ {fxxy ± 2Mo/xj/ ± fyt)fx + fxyfxx ± ftfxy\f + '^fxfyifxx ± ^o/x ± ft) = 0. (4) 

To solve the trilinear equation one has to use some prior ansatz. The variable separation 
solution (I) is resulted from the ansatz / = + ciiP + 0.2(1 + c^sM- 

For the DLWE (^, there are two sets of infinitely many symmetries and every symmetry 
possesses an arbitrary function of t or y[|l4|- That means infinitely many arbitrary functions 
of y and t can be entered into the solutions of @. So it is possible to extend the solution (H) 
to a more general one with more arbitrary functions. After finishing the detailed calculations, 
we find that the following variable separation ansatz, 

N 

f = qo + ^PiQh (5) 

i=l 

where {g^, i = 0, 1, 2, A^}, and {pi, i = I, 2, A^} are functions of {y, t} and {x, t} 
respectively, solves the trilinear equation (^) under the conditions 

N N 

lit = + qiCj)qj, Pit = (coo - uodx - dl)pi - CQi + ^(cjoPi - Cji)Pj, (6) 

i=o j=i 

for all i = 0, 1, 2, A^ and {q^, z, j = 0, 1, 2, A^, Cj, j = 1, 2, A^} being 
arbitrary functions of t. Obviously, the general ansatz (|]) will return back to the known one 
for A^ = 2, qo = ao, qi = ai, q2 = ^2?, qs = a^q, pi = ps = p, P2 = 1- 
The corresponding solution for the field v = —1] — 1 now reads 



qo + Er=i Piqi (qo + Eii p*?* 



Ue (7) 



while the quantity u is given hj u = ±{2J2iLiPixqi)/{qo + J^iLiP^qi) + ''^o- It is clear that 
in addition to a (l+l)-dimensional arbitrary function of {x, t} (one of uq and pi), {N + 
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1)(A^ + 2) — 1 arbitrary functions of t, Cij, Cj have been included in the general solution. 
Furthermore, various arbitrary functions of y and {x, t} will also be included in (|^ after 
the coupled systems of (|]) being solved. Because of the complexity, we have to leave these 
problems for our future studies. Here we just write down the simplest nontrivial case for the 
later uses. It is quite trivial that, when we fix = 1, Cij = = 0, pi = p and qo ^ aQ + go, 
the formula (|^) is simplified to 

2px {qiqoy - {ao + qo)qiy) 

V = (Ol 

{ao + qo+pqiY 

with go and qi being arbitrary functions of y and p being an arbitrary function of {x, t}. We 
can prove that the simplified quantity expressed by the right hand side of (§) does work for 
all the known models that allow the universal formula (|l]). It is interesting also that though 
we have not yet proven the validity of (0) for all the models mentioned in 0, we do have 
proven that the extended form (|^) is really valid at least for some of them such as the BKK 
system, the (2+l)-dimensional Burgers system and the (l+l)-component AKNS system. 
For instance, the quantities g = —2G = Ue and H = (In/)^ + ^0/2 with the same /, pi 
and qi determined by (|) and (||) solve the BKK system {Hty — H^xy + '2{HHx)y + '^^G^x = 
0, Gt + Gxx + 2{HG)x = Q]. 

As pointed out in [^], because of the arbitrariness of the functions appeared in the universal 
formula and its extended forms, various class of localized excitations can be constructed. 
Now we concentrate on how to find some types of folded solitary waves and foldons from the 
universal formula and its extended forms. First of all, we write a localized function, px, in 
the form 

M M 

VX = Y,fM + '^ji)^ ^ = e + $^^7,(e + ^^,i), (9) 
i=i i=i 

where f 1 < f 2 < ■ ■ -vm are all arbitrary constants and {/j, gj}^ Vj are all localized functions 

with the properties /,(±oo) = 0, (^^(±00) = Gf = consts. From the second equation of @, 

we know that ^ may be a multi-valued function in some possible regions of x by selecting 

the functions gj suitably. So the function px may be a multi-valued function of x in these 

regions though it is a single valued function of ^. It is also clear that px is an interaction 

solution of M localized excitations because of the property ^Ix^oo 00. Actually, most of 

the known multi-loop solutions are the special cases of (^. Now if we take all the arbitrary 

functions appeared in the universal formula (P and/or its slightly general one (H) (and/or 
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even its more general extended one (0)) possess the forms similar to then we can get 
various types of folded solitary waves and/or foldons. 

In Fig. 1, four typical special folded solitary waves are plotted for the quantity v shown 
by (|^) with = px, p = p^x^d^ and the functions go and qi being given in a similar way 

qsy = y = + R{0), q5= qsyVede, {6 = 0, 1). (10) 

i=i 

In (|TUp, Qsj{0), Vj and R{9) are localized functions of 9. The more detailed function selec- 
tions of the figures are directly given in the figure caption. 

To construct the possible localized excitations with completely elastic interaction behav- 
iors, especially, foldons, we should discuss the asymptotic behaviors of (§) and (||) when 
t — > ±oo. Assuming further that Fj{ujj) = Fj{^ + vjt) = J fjdx\i^.^±oo — > F^, we can 
gaze at the jth excitation to see the interaction properties among the localized excitations 
expressed by (|^) with and (|T0|). In other words, we can take ujj invariant and then take 
t ^ oo because go and gi have been fixed as t- independent. The results read 

vU 2/,(..,)[gigo,-(ao + go)g,] ^ . , + a ■(£ + v -t) (11) 

where fij = Y.i>j + Ei<j and VJ = Y.i>j Gf + Y.i<j Gf. From the asymptotic 
results (P^TD, we know that, (i) the jth localized excitation given by with (^) and (^) is a 
travelling wave moving in the velocity Vj along the negative {vj > 0) or positive {vj < 0) x- 
direction; (ii) the multi- valued properties (i.e., the structures) of the jth localized excitation 
is only determined by gj of (^; (iii) the shape of the jth excitation will be changed if 
Qj 7^ Qj , however, it will preserve its shape if Qj = Qj; (iv) the total phase shift for the 
jth excitation is Tj — Tj. 

From the above discussions, it is seen that to construct some (2+l)-dimensional localized 
soliton type excitations with completely elastic interaction behaviors becomes an easy task. 
The only thing one has to do is to select suitable localized functions Fj (and then fj) with 
Qj = Qj in (P)! The localized soliton like solutions can still possess quite rich structures. 
Especially, various types of known single-valued localized excitations are also in the catalog 
due to the arbitrariness of the functions included in (^ and (|10|). For instance, if we select 
fi and Qi of (^ and (jT^) possesses the peakon shape solutions[^ ^ and/or the compacton 
shape solutions |TH] with the property Qf = , then we can find some (2+l)-dimensional 



peakon type and/or compacton type solutions with completely elastic interaction properties. 
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Figure 1 . Four typical folded solitary wave solutions for the field v expressed 
by (H) at t = with @, (|^) and the related concrete selections are: (a) Px = 
sech^(^— t>it), qoy = sech^i], qi = 1, x = ^— 1.15 tanh(^— t>it), y = rj+lotanhr], 
Zq = 1 and ao = 8; (b) same as (a) but with Iq = —1.15 and = 4; (c) and 
Qoy are same as (a), however, qiy = secfir], x = ^ + 2 tanh(,^ — vit) + tanh^(^ — 
Vit) — kQ tanh^(^ — Vit), y = r] + 2 tanhi] + tanh^ V ~ f^o tanh'^ i], = —5.4 and 
an = 50: (d) same as (c) but with kn = 10 and an = 250. 



Figs. 2a-2c are plotted to show the possible existence of foldons which are given by @ 
with gi = 1, and (|10|) for 6 = and the concrete function selections are also given in 
the figure caption. From Figs. 2a and 2c, we can see that the interaction of two foldons is 
completely elastic. Because one of the velocities of foldons has been selected as zero, it can 
also be seen that there are phase shifts for two foldons. Especially, before the interaction, 
the static foldon (the large one) is located at x = —1.5 and after the interaction, the large 
foldon is shifted to x = 1.5. 

In summary, starting from the "universal" formula and/or its extended form, various kinds 
of localized excitations with and without completely elastic interaction behaviors can be 
constructed easily by selecting the arbitrary functions suitably according to the asymptotic 



results iuA]) . Especially, a new kinds of localized excitations, folded solitary waves and 
foldons, are investigated both analytically and graphically. The foldons may be folded quite 
freely and complicatedly and then possess quite rich structures and interaction behaviors. 
The explicit phase shifts for all the localized excitations given by (j^) with (j^) and (|T0|) have 
been given. 

On the one hand, there are a large number of complicated "folded" and/or the multi- 
valued phenomena in the real natural world. On the other hand, there is no good analytical 
way to treat these kinds of complicated phenomena. This letter is only a beginning attempt, 
to find some types of possible stable multi-valued localized excitations, folded solitary waves 
and foldons, for some real physical models including some quite "universal" systems such as 
the well known DS equation and the DLWE. 

Because the formula is valid for various (2+l)-dimensional interesting models, which 
are divergently applied in many physical fields, we believe that the foldons may be useful 
in the future studies on the complicated "folded" natural world. The more about both the 
(extended) "universal" formula and the general (or special) foldons especially its real possible 
applications should be studied further. 

The authors are in debt to thank the helpful discussions with Prof. C. Rogers, Dr. W. 
Schief, Prof. C-z Qu and Dr. H-c Hu. The work was supported by the National Outstanding 
Youth Foundation of China (No. 19925522), the National Natural Science Foundation of 
China (No. 90203001) and the Research Fund for the Doctoral Program of Higher Education 
of China (Grant. No. 2000024832) and ARC Large Grant A00103139 of the University of 
New South Wales University, Australia. 
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Fig.2a 




Fig.2c 




Figure 2. Evolution plots of two foldons for the field v expressed by (H) with 
= O.Ssech^^ + 0.5sech2(^ - 0.25t), x = ^ - l.Stanh^ - 1.5tanh(^ - 0.25t), 
qi = 1, go = sech^^, y = 9 — 2tanh6' and ao = 20 at times, (a) t = —18, (b) 
t = 7.2 and (c) t = 18 respectively. 
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